Topological invariants of time-reversal-invariant band structures 
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The topological invariants of a time-reversal-invariant band structure in two dimensions are mul- 
tiple copies of the Z2 invariant found by Kane and Mele. Such invariants protect the topological 
insulator and give rise to a spin Hall effect carried by edge states. Each pair of bands related by 
time reversal is described by a single Z2 invariant, up to one less than half the dimension of the 
Bloch Hamiltonians. In three dimensions, there are four such invariants per band. The Z2 invariants 
of a crystal determine the transitions between ordinary and topological insulators as its bands are 
occupied by electrons. We derive these invariants using maps from the Brillouin zone to the space 
of Bloch Hamiltonians and clarify the connections between Z2 invariants, the integer invariants that 
underlie the integer quantum Hall effect, and previous invariants of T- invariant Fermi systems. 
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In a remarkable pair of papers, Kane and Mele 0, Q 
proposed a Z 2 topological invariant of time-reversal- 
invariant insulators in two dimensions, showed that the 
nontrivial "topological insulator" phase created by spin- 
orbit coupling has an intrinsic spin Hall effect distinct 
from earlier proposals 0, Q , and argued that graphene 
is a viable system in which to observe this effect. (Here 
Z 2 = Z/2Z is the cyclic group of two elements.) A di- 
rect and experimentally relevant characterization of the 
invariant was given in terms of edge states at the bound- 
ary of a 2D insulator: the topological insulator has an 
odd number of Kramers pairs of edge modes, while the 
ordinary insulator has an even number. Two explana- 
tions for the invariant as a property of the bulk band 
structure with spin-orbit coupling were also given. 

We first show that time-reversal-invariant (T- 
invariant) 2D insulators have multiple Z 2 invariants 
that are directly analogous to the band TKNN integers 
or Chern numbers in the integer quantum Hall effect 
(IQHE) and can be understood using similar methods. 
This clarifies and expands upon the two original descrip- 
tions of a single bulk Z 2 invariant, which do not follow 
the standard homotopy paradigm of most topological in- 
variants in condensed matter physics, and whose connec- 
tion to the IQHE is opaque: for the case of an occupied 
pair of bands relevant to graphene, the Z2 invariant was 
explained either by invoking K-theory (recently used to 
classify Fermi surfaces pj) or by counting zeroes or zero 
regions of a certain matrix Pfaffian. 

The Z 2 invariant is first rederived for the simplest case 
(two occupied bands related by time-reversal symmetry) 
using basic homotopy theory to establish its connection 
to the usual Chern number or TKNN integer [§J of T- 
invariant systems. There is a connection between its 
absence with inversion symmetry and previous work on 
topological invariants of T-invariant Fermi systems 0. 
We then obtain results for multiple occupied bands: there 



are multiple Z 2 invariants in such systems, and perhaps 
most interestingly, four invariants per band in 3D. In 
2D, the significance of these invariants in physical sys- 
tems is straightforward: just as two IQHE states with 
the same sum of Chern numbers for occupied bands are 
adiabatically connected, two T-invariant band insulators 
are adiabatically connected if and only if they have the 
same Z 2 sum of individual Z 2 invariants. 

The approach in this paper can be outlined as follows. 
The basic objects of homotopy theory are the homotopy 
groups 7r„(M) that describe equivalence classes under 
smooth deformations of mappings from the sphere S n 
to a manifold M. A band structure can be thought of 
as a map from the Brillouin zone (a torus rather than 
a sphere) to the space of Bloch Hamiltonians, assuming 
that the Hilbert space is the same for all points in the 
Brillouin zone. The T symmetry means that the effective 
Brillouin zone (EBZ), a set of points for which the Bloch 
Hamiltonians can be specified independently, is nearly a 
sphere: we consider "contractions" that extend a map- 
ping from the EBZ to one from the sphere. 

A given mapping from the EBZ has an infinite set of 
topologically inequivalent contractions in this sense. We 
consider a mapping from the EBZ together with the set 
of all of its contractions, which is invariant under smooth 
deformations of the original mapping. There are then 
exactly two equivalence classes, for the case of two occu- 
pied bands connected by T, which correspond to the ordi- 
nary insulator and the topological insulator. In this sim- 
plest case, an ordinary (topological) insulator becomes 
the equivalence class of all mappings from the sphere with 
even (odd) Chern numbers. This construction requires 
only standard homotopy results, makes no assumptions 
about the details of the band structure or the existence of 
additional commuting operators such as spin, and gener- 
alizes directly to multiple bands and higher dimensions. 

For 2D T-breaking systems, in a nondegenerate band 
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structure each band is associated with a TKNN integer 
that is invariant under smooth perturbations of the Bloch 
Hamiltonians. More precisely, if the Hilbert space of the 
Bloch Hamiltonians is finite and of dimension n, there are 
n — 1 independent integer-valued invariants |8| because 
the invariants sum to zero. There are a number of ways 
to obtain the TKNN integers, e.g., as an integral over the 
magnetic Brillouin zone (a torus) of the Berry flux, using 
the explicit wavefunction ipi for band i, or as an integral 
using the projection operator || 
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Tr(dppdp), Pj = |^X^I- (1) 



BZ 



Here dPi — dxd x Pi + dyd y Pi and dxdy = —dydx. A 
powerful way to understand these integer invariants || 
is by considering mappings of the torus to Bloch Hamil- 
tonians, which arc assumed to be nondegenerate and act 
upon the same Hilbert space. Denoting the space of such 
Hamiltonians as A4, the existence of the TKNN integers 
follows from the first two homotopy groups 



TTl (M) = 0,7T 2 (M) = Z" 



(2) 



where the second formula indicates n — 1 copies of the in- 
finite cyclic group Z. The second result follows [8| from 
the same exact sequence as used in the theory of topo- 
logical defects 9]. For a pair of bands i, j that are possi- 
bly degenerate with each other but with no other bands, 
there is a single integer-valued invariant || obtained by 
replacing Pi in with Pjj = Pj + Pj ■ 
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FIG. 1: The topology of the effective Brillouin zone (EBZ): 
if the original Brillouin zone is the torus in (a), then T- 
invariance reduces the independent degrees of freedom to live 
on the manifold in (b). Points on the boundary circles that 
are connected by horizontal lines are conjugate under T; the 
points r, A, B, C are self-conjugate, and their Bloch Hamilto- 
nians are therefore in the even subspace Q. 

Now consider consequences of invariance under the 
time-reversal operator T. For fermions, T 2 = — 1 and 
T is represented by an antiunitary operator in the 
Hilbert space of Bloch Hamiltonians. Time-reversal con- 
nects both pairs of points in the Brillouin zone (fc, — fc) 
and the associated Bloch Hamiltonians: 



Fig. 1 shows the original toroidal Brillouin zone and the 
EBZ: specifying the Hamiltonians for all points on the 
EBZ determines them everywhere. The Bloch Hamilto- 
nians can be specified independently except at the bound- 
aries, where points are related by T. Clearly points at 
which k = — fc, such as T, A, B, C in Fig. 1, are special: 
at these points the Bloch Hamiltonian commutes with 
O. This "even" subspace in the language of Ref. Q has a 
natural interpretation as the set of symplectic or quater- 
nionic Hamiltonians Q: we denote this subspace, with 
the additional assumption of no degeneracies other the 
two-fold degeneracies required by T, as Q. T-invariance 
requires an even number of bands In. 

In general, a T-invariant system need not have Bloch 
Hamiltonians in Q except at these special points as long 
as inversion symmetry is broken, so that H(k) ^ H(—k). 
At first glance, there is no obvious topological invari- 
ant for a degenerate band with time-reversal symmetry, 
because the Chern number for the whole Brillouin zone 
vanishes. It is simplest to see this for the projection op- 
erator Pi corresponding to a single nondegenerate band: 
using an explicit representation Pj = |?/>j}(?/>i|, 

rii = ±- [ Tr(dPiPidPi) 

M JBZ 

(d x ipi\d y ^i) + (d x <di/ji\d y Qipi) 
-Im / {d x ipi\d y ipi) + (d v 4>i\d x ipi) =0.(4) 
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H(-k) = QH{k)Q- 



(3) 
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A similar argument applies for two possibly degener- 
ate bands, but there is nevertheless a topological invari- 
ant 0], and in general multiple topological invariants. 

As a quick example of homotopy arguments, suppose 
that inversion symmetry is unbroken, which implies that 
the Bloch Hamiltonians are everywhere in Q. Any map- 
ping of the torus T 2 to Q with the condition that points 
k and —k go to the same point is determined by its be- 
havior on one point from each (fc, — fc) pair, and the EBZ 
under this condition is topologically identical to a sphere 
(stitching T-conjugate points together in Fig. lb): the 
classes of such mappings are given by 
Hence there is no homotopy invariant for 2D band struc- 
tures with both T and inversion symmetry, although 
^i(Q) 7^ and higher-dimensional invariants can exist. 

Now consider possible invariants without inversion 
symmetry. We seek to classify mappings from T 2 to the 
space of Hamiltonians C that have at most two-fold de- 
generacies within a pair of T-related bands and behave 
under O as specified above. Such a mapping is deter- 
mined by a mapping from the EBZ, i.e., a mapping from 
the cylinder C to C with certain conditions on the two 
circular boundaries reflecting time-reversal. The image of 
a boundary point must be the Q conjugate of the image 
of the point on the same boundary related by k <-» — k. 

If the same Hamiltonian occurred at all points on the 
boundary, then the topology of the cylinder becomes that 
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of the sphere, and the degenerate Chern numbers |8| are 
integer-valued invariants for mappings from the sphere 
to C that give one integer for each possibly degenerate 
pair since 7r 2 (C) = Z™ -1 . For the simplest case of two 
occupied possibly degenerate bands, the invariant com- 
bination reduces to n\+ri2 if the bands are nondegenerate 
and hence have separate Chern numbers n\ and rt 2 . 

We show first for the simplest case of one pair of bands 
that any mapping from C to C, even if the elements at a 
boundary are not all the same, can be smoothly deformed 
("contracted") to one in which the boundary elements 
are identical to an arbitrary reference element Qq 6 Q 
(Fig. 2a); the resulting map from the sphere has a well- 
defined Chern number. It is required that at each stage of 
the contraction, the boundary has the same conjugacy of 
points under T as in the original boundary. This guaran- 
tees that two maps from the EBZ that can be contracted 
to maps from the sphere with the same Chern number 
are homotopic (deformable to each other). 

Then it is shown that different contractions differ by an 
arbitrary even Chern number, so that there are only two 
invariant classes according to whether the Chern num- 
ber with any contraction is an odd or even integer. The 
existence of at least one contraction follows immediately 
from 7Ti (C) = 0: contract one side of the boundary circle 
on the point Qo, keeping the top and bottom points in Q, 
and then determine the other side's behavior by T con- 
jugacy. The fact that different contractions differ by an 
even Chern number is less trivial and gives an unexpected 
picture of how the Z 2 invariant arises: as an embedding 
of one Z homotopy group in another. Whether a band 
is "odd" or "even" is directly computable by using any 
contraction to define the Chern number integral. 

A contraction in the sense above gives a smooth map- 
ping from a cylinder to C, f(6,X), < A < 1 such 
that for each A, f(9, A) and f(2n — 9, A) are T conju- 
gates (which implies that /(0, A) and f(ir, A) are in Q). 
Now f(9, 0) should agree with the initial specification of 
Bloch Hamiltonians, while f(9, 1) = Qo is constant. If 
both boundaries are contracted, the resulting sphere has 
a well-defined Chern number for each pair of bands. 

In fact, many topologically inequivalent contractions 
exist, and these different contractions are responsible for 
reducing the integer-valued invariants on the sphere to 
Z 2 invariants on the EBZ. Let /1 and / 2 be two different 
contractions. Then define a mapping g(9,X), which by 
composition changes from contraction fi to / 2 (Fig. 2b): 



9(0, X) 



/i(0,l-2A) if < A < 1/2 
/ 2 (0,2A - 1) if 1/2 < A < 1 



(5) 



Although the domain of the mapping g is topologically 
a sphere because the circles at both ends go to the same 
point, g differs in its T symmetry from the contracted 
half of the Brillouin zone, which in its interior has no T 
symmetry relating different values of 9. Let the cylinder 



in the definition of g have coordinates 8 £ [0, 2tt), A £ 
[0, 1]; then points (8, A) and (2ir — 9, A) are T-conjugate. 
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FIG. 2: (a) Contracting the extended Brillouin zone to a 
sphere, (b) Two contractions can be combined according to 
© to give a mapping from the sphere, but this sphere has a 
special property: points in the northern hemisphere are con- 
jugate under T to those in the south, in such a way that 
overall the Chern number must be even. 

It follows that g has arbitrary even Chern numbers, 
and the change in the Chern numbers of Ub induced by 
changing from contraction fx to contraction / 2 is 



An l B = 2n l s , 



(6) 



where i indexes band pairs. This can be verified in two 
steps: map the equator of the sphere S to the constant 
element E , which is possible since ni(Q) = and topo- 
logically unique since 7r 2 (Q) = 0, then note that each 
hemisphere has a well-defined Chern number and that the 
Chern numbers of the two hemispheres are equal, rather 
than opposite as in the case of the original Brillouin zone. 
The reason for this equality can be understood easily in 
the cylindrical coordinates above, where the equator is 
at 9 = it and 6 — 0. The identification under T of 9 and 
2n — 9 means that d9 changes sign between a point and 
its time-reversal conjugate, but dX does not, giving an 
additional change of sign in equation Q . 

The same topological argument applies to n pairs of 
bands. Since there is one integer invariant for each such 
pair with a zero sum rule, there is one Z 2 invariant for 
each pair, with an even number of "odd" band pairs. 

We now give the generalization to three-dimensional 
Brillouin zones, where there are significant differences be- 
tween Z 2 invariants and the 3D integer-valued TKNN in- 
variants 0, : there are four independent Z 2 invariants 
per pair of bands, even though there are only 3 Chern 
numbers for a pair of degenerate bands. An alternate 
way to obtain the 2D result, which is more useful for 3D, 
is by considering contractions to a torus, rather than a 
sphere. The set of mappings from the Brillouin zone T 3 
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to C is determined by three ordinary Chern numbers since 
7T3(C) = 0: the three integers correspond to the xy, yz, 
and xz planes The mappings from any two parallel 
planes in T 3 to C have the same Chern number. 

Now consider possible Z2 invariants for a time-revcrsal- 
invariant system in 3D. Suppose that the Brillouin zone 
torus is x,y,z e [—1, 1], and construct an EBZ by taking 
the part of the torus with z > 0. Then slices at constant 
z for < z < 1 have the topology of the torus T 2 , while 
at z = and z = 1 there are additional T constraints 
that reduce the degrees of freedom to the 2D BZ. The 
T constraint means that the xy Chern number is zero in 
every slice, while the xz and yz integers are arbitrary. 

The EBZ boundaries at z = and z = 1 are charac- 
terized by one Z2 invariant each: one boundary may be 
even while the other is odd because the boundary slices 
have the same Chern number (zero) and thus are homo- 
topic as maps to C in the EBZ interior. Once the two Z2 
invariants at the boundaries are fixed, two contractions 
of an orig inal 3D EBZ to the torus T 3 differ by two even 
Chern numbers, one for the xz slices and one for the yz 
slices. These are even because a slice of the contraction 
has the same symmetries that force even Chern numbers 
in the 2D case. There are four additive Z2 invariants per 
band pair in 3D and 16 insulating phases. 

Going back to the torus T 3 , there are 6 inequivalent 
planes that have the symmetries of the 2D BZ and hence 
have a Z2 invariant. Only four of these are independent, 
consistent with the counting above: if the six invariants 
with values ±1 are xq, x\, yo, yi, zq, z±, for the planes x = 
and x = ±1 and so on, then there are two relations 

xqXi = y yi = z z\. (7) 

For example, a model can be designed on the 3D sodium 
chloride lattice to reduce to a previously introduced 
2D square lattice topological insulator in the k y = or 
k z = planes: it has a phase with yo = y± = z$ = Z\ = 
— 1, xo — x\ = 1. In 2D, both insulating phases can 
be realized in models where there is a conserved quan- 
tity (e.g., S z ) that allows definition of ordinary Chern 
integers. A new feature in 3D is that the phases with 
xqXi = z/oj/i = z a z i = — 1 cannot be realized in this way. 

Since the results here are for Hilbert spaces of ar- 
bitrary dimension, they can be applied to many-body 
problems with an odd number of fermions [7j if there 
are two periodic parameters in the Hamiltonian that are 
connected by time-reversal in the same way as the mo- 
mentum components (k x ,k y ). Our derivation of Z2 in- 
variants helps explain spin Hall edge states Q]: in 2D, 
just as Chern number predicts the number of edge states 
in the IQHE [llj, the class of even (odd) Chern num- 
bers in the bulk corresponds to edges with an even (odd) 
number of Kramers pairs of modes. Note that transport 
by these edge states will receive power-law, rather than 
exponential, corrections in voltage or temperature. 



The existence of a 2D Z2 invariant was also obtained by 
Haldane . The bulk-edge connection has been derived 
when ordinary Chern integers are defined BQEl (see 
also Ref. Hit) . Recently preprints appeared on the 2D 
Z2 invariant as an obstruction [Tij or as a noninvariant 
Chern integral plus a formal integral on the EBZ bound- 
ary that is defined up to addition of an even integer |l6| . 
Defining this integral is equivalent to our prescription of 
choosing any contraction to define a Chern integer. To 
our knowledge, the full counting of Z2 invariants in 2D 
or 3D has not been obtained before. 

Although whether graphene has sufficiently strong 
spin-orbit coupling to realize the topological insulator 
has been debated [ItI Hj| . it is hoped that the under- 
standing of Z 2 invariants developed here will stimulate 
searches in a wider class of materials, especially in 3D. 
The understanding of bulk Z2 invariants in this paper, 
combined with the stability of total Z2 to interactions 
and scattering at the edge 0, H3] , shows that the topo- 
logical insulator is a robust phase of matter with a deep 
connection to the quantum Hall effect. 
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